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Abstract 

A new development of the "monodromy transform" method for analy- 
sis of hyperbolic as well as elliptic integrable reductions of Einstein equa- 
tions is presented. Compatibility conditions for some alternative repre- 
sentations of the fundamental solutions of associated linear systems with 
spectral parameter in terms of a pair of dressing ("scattering") matrices 
give rise to a new set of linear (quasi-Fredholm) integral equations equiv- 
alent to the symmetry reduced Einstein equations. Unlike previously de- 
rived singular integral equations constructed with the use of conserved 
(nonevolving) monodromy data on the spectral plane for the fundamen- 
tal solutions of associated linear systems, the scalar kernels of the new 
equations include another kind of functional parameters - the evolving 
("dynamical") monodromy data for the scattering matrices. For hyper- 
bolic reductions, in the context of characteristic initial value problem these 
data are determined completely by the characteristic initial data for the 
fields. In terms of solutions of the new integral equations the field com- 
ponents are expressed in quadratures. 



1 Introduction 



In General Relativity in a number of physically significant cases the dynamical 
part of Einstein equations, being restricted to space-time geometries admitting 
two-dimensional Abelian isometry groupfj reduces to the nonlinear integrable 
systems. Among these are the Einstein equations for vacuum gravitational fields 
[p|— IpH, the Einstein equations for space-times with a stiff matter fluid fl^] , the 
clcctrovacuum Einstein - Maxwell field equations |l^, [n], 13 1, the Einstein - 
Maxwell - Weyl equations for gravitational, electromagnetic and massless spinor 
fields Qj, as well as some string theory induced gravity models, e.g., the Einstein 
- Maxwell equations with axion and dilaton fields Q| - . Accordingly to the 
type of the contemplated two-dimensional space-time symmetry (determined 
by the signature of the metric on the orbit space) the reduced equations can 
be either of the hyperbolic or of the elliptic types. By now a theory of these 
equations has been developed and discussed in many aspects (for details and 
references see, for example, |l6|, |lj). 

The so called monodromy transform jl6|, [l8|, [llj provides some general and 
fairly simple base for a description in a unified manner of all mentioned above 
integrable reductions of Einstein equations. Similarly to the well known inverse 
scattering method (the "scattering transform" ) , this approach begins with a 
representation of the dynamical part of reduced Einstein equations as the in- 
tegrability conditions of an overdetermined linear system of a special structure 
containing a spectral parameter. The analysis of the constructed linear systems 
showed a common important property of the evolution of fields described by all 
mentioned above integrable reductions of Einstein equations. This property is a 
conservation of the monodromy structure on the spectral plane of the normalized 
fundamental solutions of associated linear systems. This monodromy structure 
is defined here as a set of linear transformations which relate this fundamental 
solution with its analytical continuations along the closed paths surrounding 
each of the singular points of this solution on the spectral plane. The matrices 
of these linear transformations possess some special algebraic structure with a 
small set of independent components - the functions of the spectral parame- 
ter which constitute what we call here as monodromy data.0 A remarkable 
properties of the defined monodromy data are that they (i) are functions of the 
spectral parameter only, (ii) exist for any analytical local solution of reduced 

1 The most known and elegant form of the reduced equations is written in terms of complex 
potentials — the Ernst potentials. These are the vacuum Ernst equation for one complex 
potential or electrovacuum Ernst equations, i.e. a coupled system of two similar equations 
for two complex Ernst potentials fcj] which have been derived by F.J.Ernst originally for 
stationary axisymmetric fields. For the hyperbolic case these equations have very similar 
forms {a. A generalized form of these equations arises, for example, in the presence of Weyl 
spinor field M. 

2 For example, for any vacuum gravitational field with the supposed space-time symmetry 
the monodromy data consist of two functions of the spectral parameter, while for elctrovacuum 
fields we have four such functions. 
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Einstein equations and (iii) characterize uniquely every analytical local solu- 
tion of the field equations. All these means, that these monodromy data can 
be considered as a new set of the field variables for which the field equations 
become trivial (being the conserved quantities, these monodromy data do not 
possess any evolution) or as new "coordinates" (instead of usual field variables) 
in the infinite dimensional spaces of local solutions of the integrable reductions 
of Einstein equations. Just using of such "coordinate transformation" explains 
why this approach was called as the monodromy transform. 

In this way the "direct" and "inverse" problems of the monodromy transform 
suggest themselves naturally. For solution of the first of them, i.e. for calcula- 
tion of the monodromy data for a given solution of reduced Einstein equations, 
we have to find a fundamental solution of the associated linear system with the 
coefficients corresponding to given field components and to determine than its 
monodromy data on the spectral plane. The solution of the inverse problem, i.e. 
the construction of a solution of Einstein equations for given monodromy data 
functions, reduces to the solution of a special integral equation form of the field 
equations derived in H, [l9| . These are the linear singular integral equations 
whose scalar kernels and right hand sides are expressed in terms of the men- 
tioned above nonevolving monodromy data. (Equivalent regularizations of these 



equations have been derived in MA 16, Gl|.) These integral equations always 



possess a unique solution for any given monodromy data functions holomorphic 
in some local regions of the spectral plane. The corresponding local solution 
of Einstein equations can be calculated in quadratures whose integrands are 
expressed in terms of the solution of the mentioned linear integral equations. 

This approach suggests different applications. The first of them is a direct 
construction of local solutions of Einstein equations. In this case the actual 
problem is to chose some monodromy data which allow to solve explicitly the 
integral equations with the corresponding kernels and then to calculate the cor- 
responding quadratures for the field components. Fortunately, for very large 
classes of the monodromy data all these can be realized, and infinite hierarchies 
of families of exact solutions with any finitely large number of free parameters 
can be calculated in elementary functions. The examples of such hierarchies 
arc the solutions with arbitrary analytically adjusted rational monodromy data 
23| and the solutions for colliding plane waves and inhomogeneous cos- 
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mological models with analytically not adjusted but also rational monodromy 
[ 24j . Each of these hierarchies extends considerably the hierarchies of vac- 
and electrovacuum ||l3f multisoliton solutions provided the Minkowski 
space-time is chosen as the background for solitons. 

Another application of the same approach is a formulation of general schemes 
for solution of some initial or boundary value problems for the gravitational 
and some matter fields with two-dimensional space-time symmetries. The basic 
idea is that any given initial data for a characteristic initial value problem or 
a Cauchy problem for integrable hyperbolic reductions of Einstein equations 
as well as a boundary data for some boundary problems for an elliptic (e.g., 
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stationary axisymmctric) integrable reduction of Einstein equations allow to 
calculate the corresponding monodromy data. The conserved (i.e. nonevolving) 
character of these data allows to identify them with the monodromy data for the 
sought-for solution of the initial or boundary value problem under consideration. 

In the present paper a derivation of a new form of integral equations equiva- 
lent to each of the mentioned above integrable reductions of Einstein equations 
is presented. For this we introduce a specific alternative representations of 
the fundamental solutions of the associated linear systems in terms of pairs of 
"scattering" matrices dressing some partial values of these fundamental solu- 
tions ("in-states"), each depending on one of two coordinates on the orbit space 
and on the spectral parameter. These scattering matrices were found to possess 
a specific analytical structure on the spectral plane, where each of them is char- 
acterized by two (algebraic in the absence of a Weyl spinor field) branchpoints 
and finite jump on the cut which joins these points. (It is useful to recall here 
that the fundamental solution of any of the associated linear systems under 
consideration in general possess four such branchpoints which are joined in our 
construction by two nonintersecting local cuts.) The consistency conditions for 
these alternative representations of solutions of associated linear systems give 
rise to linear (quasi-Frcdholm) integral equations interrelating some fragments 
of the algebraic structures of the dressing matrices on the cuts. The scalar 
kernels of these equations and their right hand sides include functional param- 
eters which characterize the monodromy properties of the dressing (scattering) 
matrices mentioned above. Unlike the nonevolving monodromy data for funda- 
mental solutions of associated linear systems, these ("dynamical") monodromy 
data evolve and their evolution is determined by a coordinate dependence of 
the "in-states". In terms of solutions of these new integral equations all of 
the field components and the Ernst potentials, characterizing the solutions, are 
determined in quadratures. For hyperbolic reductions, in the context of charac- 
teristic initial value problem these "in-states" can be identified easily with the 
characteristic initial data for fundamental solution of associated linear system 
which are determined completely by the characteristic initial data for the fields. 

For simplicity, we restrict all considerations following below by hyperbolic 
as well as elliptic reductions of vacuum Einstein equations and electrovacuum 
Einstein - Maxwell equations only, because there are no any principal difficulties 
for the realization of the constructions suggested below for all other mentioned 
above integrable reductions of Einstein equations. 

2 Reduction of Einstein equations 

In this section we describe basic definitions used below and present in unified 
notations the dynamical parts of reduced vacuum Einstein equations and elec- 
trovacuum Einstein - Maxwell equations in one of their most compact forms — 
in the form of the Ernst equations. 
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2.1 Metric and electromagnetic potential 

The components of a 4-dimensional space-time metric and a 1-form J> of com- 
plex electromagnetic potential for a self-dual Maxwell 2-form in a space-time 
admitting 2-dimcnsional Abelian isometry group can be considered in the form 

ds 2 = g^dx^dx" + g ab dx a dx b , $ = $ a dx a (1) 

where /x, u, . . . = 1, 2; a, b, . . . = 3,4; g^ v , g a b and $ a depend on the coordinates 
x 1 and x 2 only. We consider metric g^ locally in a conformally flat form and 
parametrize the nonzero field components in ([l]) by scalar functions 

(2) 

where / > 0, a > 0, H > 0; e, eo, ei, £2 = ±1 are the sign symbols whose choice 
should provide the Lorentz signature of the metric (|l|). This condition is equiva- 
lent to the relation £\E2 — — e. The sign symbol e determines both, the signature 
of metric g a b on the 2-dimensional orbits of the space-time isometry group and 
the signature of the conformally flat metric g^ on the two-dimensional orbit 
space of this group. We shall refer to the cases e = 1 and e = — 1 as the 
hyperbolic and elliptic ones respectively. 

The function a(x 1 ,a; 2 ) in (0) satisfies the identity det \\g a b\\ = ea 2 , which 
means that this function characterizes a measure of area on the orbits of the 
isometry group. For all known integrable reductions of Einstein equations con- 
sidered here a(x 1 , x 2 ) is a harmonic function, and this permits to define a func- 
tion (3 as its harmonic conjugation as follows 

^=o. ^ = (^ M F »=( e O 

d,0 = -ee^a, V \ e 2 ) £ » { -e 2 J 

The functions (a,/3) constitute a set of "geometrically denned" coordinates 
which can be used instead of unspecified coordinates (x 1 ,^ 2 ).^] However, we 
use instead of coordinates a and /3 their linear combinations (£,T]): 

I Z = i + i. a where j = J 1 f ° r € = l 
y i] = (3 — ja for e = — 1 

which are two real null cone coordinates in the hyperbolic case (e = 1) or 
complex conjugated to each other coordinates in the elliptic case (e = —1). 

The existing gauge freedom permits without any loss of generality to impose 
on the metric and electromagnetic functions the normalization conditions 

= 1, n(£o,r?o) = 0, $(&,»») = (3) 
where (£oj?7o) are the coordinates of a chosen "reference" point Pq. 

3 In the stationary axisymmetric case these geometrically defined coordinates are known as 
cylindrical Weyl coordinates (p,z). 
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2.2 The Ernst equations 



The dynamical part of the space-time symmetry reduced electrovacuum Einstein 
- Maxwell field equations for hyperbolic as well as for elliptic cases can be 
expressed in the form of the Ernst equations. In our notations these equations 
and the linear equation for a, can be written as 

(Re £ + $$V(<9 M + oT 1 d lll a)d v £ - rf^E + 2$d^)d v £ = 
(Re £ + + or^d^dv® - if u {d^£ + 2$<9 A1 $)<9„$ = (4) 

rj^d^a = 

where £(x 1 ,x 2 ) and ^(a; 1 ,^ 2 ) are complex Ernst potentials, and for vacuum 
$ = 0. The Ernst potential £ is defined by the expressions jj], ||: 

Re£ = e a H - 4>$, d^QmS) = -oT^e/o^fi + - 

while the electromagnetic Ernst potential $ coincides with the corresponding 
component of the self-dual electromagnetic potential shown in The Ernst 
equations (Q) are quasilinear equations of the hyperbolic type for e = 1, and they 
are of the elliptic type for e = — 1. In the coordinates £, r] we put a — (£ — 77) / 2 j 

and rf v = f ^ J ~\ with an appropriate choice of the sign of the conformal 

factor /. 



3 The monodromy transform 

In this section we recall the basic constructions of the monodromy transform 
approach developed in ]l8l |l9[ [l6| and applicable in a unified manner to the 
analysis of all integrable reductions of Einstein equations mentioned above and, 
in particular, to the Ernst equations ([|) for vacuum and electrovacuum fields. 
The basic idea of this approach is the using of a specially defined functional 
parameters (characterizing every local solution) as new "coordinates" in the in- 
finite dimensional space of local solutions of reduced Einstein equations instead 
of usual field variables. A remarkable property of these nonevolving parame- 
ters, depending on the spectral parameter only and being interpreted as the 
monodromy data of the fundamental solutions of associated linear systems, is 
that the field equations do not impose any constraints on these parameters and 
therefore, such "coordinate transformation" solves completely the field equa- 
tions. Thus, the solution of reduced Einstein equations becomes equivalent to 
solution of the inverse problem of such "coordinate transform" (called as "mon- 
odromy transform"), and solution of this problem turns out to be equivalent to 
solution of some linear singular integral equations with scalar kernels. 
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3.1 Associated linear systems 



Among various gauge equivalent linear systems with constant or coordinate de- 
pendent spectral parameters (see |^5| for more details) we chose the Kinnersley- 
like linear system whose appropriately normalized fundamental solution seems 
to possess the most simple general analytical structure on the spectral plane. A 
complete set of matrix relations, equivalent to the reduced Einstein - Maxwell 
equations, can be expressed in terms of the four unknown matrix functions 
which are 2 x 2-matrices for vacuum gravitational fields or 3 x 3-matrices for 
gravitational and electromagnetic fields (for other integrable cases see @): 



(5) 



where w is a spectral parameter and £, rj are the defined above real (the hyper- 
bolic case) or complex conjugated to each other (the elliptic case) coordinates. 

The first group of constraints imposed on the matrix functions (^J) consists 
of two systems of linear differential equations for * with algebraic constraints 
imposed on their (also unknown) matrix coefficients 



2i{w = u (£,??)* 

2i(w - ri)d v * = V(£,7y)* 



rankU = 1, 
rankV = 1, 



trU = i 
tr V = i 



(6) 



The second group of constraints provides the existence for these linear systems 
of a common Hcrmitian matrix integral of a special structure 



wjH = WoH 



dW 

dw 





' 


1 




fj = ( 




















(7) 



where "'" is a Hermitian conjugation such that \E ,t (£, 77, u>) = ^ T (£_,ri,w) and 
Wo(ro) is an arbitrary Hermitian matrix function of the spectral parameter. In 
a vacuum case the third rows and columns of all matrices should be omitted.^] 
The third group consists of pure gauge conditions imposed without any loss 
of generality on the values of * at the chosen reference point Pq and Wo(tt): 

=1. WoH =4i(w-/3o)tt + diag(-4ee ao,-4eo,l) (8) 
where a = (£ Q - r] )/2j and /3 = (£0 + 7 ?o)/2. 



3.2 Field components and potentials 

The conditions (|^) - (|^) contain all information about the specific structures 
of U, V and W which these matrices should possess to be correctly expressed 

4 In all previous author's formulations of these groups of conditions, there were included 
also an additional condition for electrovacuum case, that the lower right element of W should 
be equal to 1. However, this condition turns out to be pure gauge one, and it can be satisfied 
by an appropriate choice of the normalization conditions. 
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in terms of the components of metric and electromagnetic potential, about the 
Ernst potentials and their relations to the metric and electromagnetic potential 
components, about the Ernst equations and all properties of the function a de- 
scribed above (see [Jl9| for some explanations). So, the relations given below are 
not some additional constraints imposed on the introduced auxiliary functions, 
but they follow directly from (||) - (||). In turn, the relations, given below in 
this subsection, take place in general, for any local solution of the reduced elec- 
trovacuum Einstein - Maxwell equations. In particular, the matrices U and V 
always possess the structures 

U = T v ■ U • T7j x 



-l 





( 1 
















Tu = 


P+ 


1 


i) 




(:': 


1 


i) 





















where the scalar functions p± = fl =p , q± — 2$ — 2$p± and 

e H 




Here £ and $ have to be identified with he Ernst potentials which character- 
ize any local solution of reduced vacuum Einstein or electrovacuum Einstein - 
Maxwell equations. The matrix function W is linear with respect to the spec- 
tral parameter w, and its components are algebraically expressed in terms of 
the metric and complex electromagnetic vector potential components: 

/ _4/ 1 a6 + 4$ a$6 _ 2 $ a \ 

W = 4i(w - + G, G = _, 

\ -2$ h 1 J 

where h ab = ea 2 g ab , g ab is the matrix inverse for the metric components g a b, 
the column <I> a = ( $, — $ ) T with the superscript T meaning a transposition. 
Accordingly to (||) and (||) , Wo (w) coincides with the value of W at the reference 
point (the normalization point) Pq, namely Wq(w) = W(£o, ?7o, w). 



3.3 The monodromy structure of rj,w) 

For any local solution of the field equations under consideration the correspond- 
ing solution of (§|)-(|j) for rj, w), considered as a function of w for given ^ and 
r\ running some local domains near their initial values £o and ijo respectively, pos- 
sesses a number of universal analytical properties on the spectral plane Jl^, [l6| . 
In particular, the components of T],w) and 1 4 r_1 (^, r], w) are holomorphic 
everywhere on the spectral plane (including w = oo, where r),w = oo) = I) 
outside two nonintersecting local cuts L + and L_ which endpoints are two 
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ordered pairs of the branchpoints (w = £o, w = £) and (w = t]q, w = if) respec- 
tively. In the hyperbolic case these cuts are represented by two nonoverlapping 
segments of the real axis on the w-plane, while in the elliptic case these cuts are 
located symmetrically to each other with respect to this axis in the upper and 
lower half-planes respectively plj. This description of analytical properties of 
iff(£,r],w) we conclude presenting general expressions for the local structure of 
^ and near the branchpoints and the cuts L± which join them |18L pj|: 



La 



*(£, r), w) = A;> r), w) ® k+(w) + M+(i, t], w) 
*- x (e, r h w) = A+l+H ® cp +(£, V , w) + N+(£, 77, w) 
77, w) = AI 1 ^ (£, 77, u;) ® k_ + M (£, rj, w) 
77, «;) = A_l_ (w) ® <p _ (£, 77, 77j) + N_ (£, jj, tu) 



A_ 











777 






- % 



(9) 

where A+ and A_ are holomorphic at w = 00, provided X±(w — 00) = 1, and 
A+ and A_ possess the jumps on L + and L_ respectively; all fragments of these 
local structures of \& and \I> _1 , i.e. each of the row and column vectors k±(iu) 
and l±(w), the row and column vectors xp ± 77, zw) and cp ±(£, 77, ttj) and the 
matrices M±(£, 77, w) and N±(£, 77, u>) are regular (holomorphic) on the cuts L + 
or L_ respectively their suffices, and the algebraic relations are satisfied on L±: 



L+,: 



k+ • N+ = 0, M+ • 1+ = 0, cp + • M+ = 0, N+ • tp + = 
k_-N_=0, M_-1_ = 0, cp_-M_ = 0, N_-tp_=0 



(10) 



The local structure (|9|) with the constraints (Jlfj) allow to clarify the monodromy 
properties of tj,w) near its branchpoints. This structure is determined by 
the monodromy matrices T + or T characterizing the linear transformations 
which relate the matrix *E r (£, Tj,w) and its analytical continuations along, say 
the clockwise directed paths t + and joining the left edge of the cut L + or 
L_ respectively with its right edge: 



*(£,ti,w) -Z+*(£,ti,w)-T±(w), T±H =1-2 



l±(w) ® k±(w) 
(l±(w)-k ± (w)) 



(11) 



For derivation of these expressions, besides (||) and (|10|), we have used the 

properties A + — A+ and A- — A_. One can observe easily that each of 
the monodromy matrices T±(w) satisfies the identity T^(uj) = I.Q 

As it was shown in fr§|, |l9| ], the linear relations between kj_(iy) and l±(w) 
implied by the conditions (f?1) allow to express the components of the monodromy 

5 We stress the point, that these properties take place for vacuum and electrovacuum cases 
only. These don't hold in the presence of the Weyl spinor field, wherj t i e lo cal structures (g) 
become more complicate and the branchpoints can be nonalgebraic 118, flq|. 
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matrices ( pi] ) in terms of the "projective" vectors k±(iu). These relations and 
the affine parametrization for k± (w) of the form 

l ± (w) = 4e (w - ^ )(w - rj Q )W 1 (w) -k^iw), k±(w) = {1, u±(w), v±(w)} 

(12) 

allow also to express all components of the monodromy matrices T± (w) in terms 
of four scalar functions u±(ui) and v±(w;) which we call as the monodromy data. 
These data are conserved, i.e. they are coordinate independent, and they can be 
determined (at least in principle) for any local solution of the reduced vacuum 
or electrovacuum field equations. For vacuum v±(w) = and therefore, any 
vacuum solution is characterized by two functions u±(w) only. 

3.4 The inverse problem of the monodromy transform 

Another important property of the nonevolving monodromy data is that they 
characterize unambiguously every local solution. For a given monodromy data 
the corresponding local solution can be calculated in quadratures in terms of 
solution of the linear singular integral equations which scalar kernels are con- 
structed using a given monodromy data. It is remarkable, that for any choice 
of the monodromy data u+(u;), v + (u>) holomorphic near the point w — £o 
and u_(u>), v_(u>) holomorphic near the point w — 770 the solution of these 
singular integral equations always exists and it is unique. These integral equa- 
tions, solving the inverse problem of such monodromy transform, are equivalent 
to the integrable reductions of Einstein equations. In |l6|, ^j] these equations 
have been presented in two alternative forms and together with their equivalent 
regularizations. Here we recall the singular forms of these equations only: 

w^o (e c)dc=k(T) ^/^rO (13) 

L L 

with scalar kernels of the following structures (the arguments £, r] are omitted): 

£(r,C) = - W C (k(r) -1(0). £(t,0 = -[A-lcWC) -l(r)) (14) 

The coordinates £ and 77 enter the integral equations ( |l3| ) as parameters which 
determine the location of the endpoints of the integration paths and as argu- 
ments of the A-multipliers in the kernels (|l4|). The expressions [A]^ in ([l4]) mean 
the jump (i.e. a half of the difference between the values of the function on the 
left and right edges of a cut) of the function A at the point ( 6 L. It seems 
useful to recall here the way for derivation of the basic equations (|l^) which 
was suggested in |1| and to explain with more details their structure. 

The described above general analytical properties of * and 'J' -1 on the 
spectral plane permit to represent these matrix functions as Cauchy integrals 
over the cut L = L + + L_ where the integrand densities are the jumps of these 
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matrix functions on L. These jumps are represented by the first terms in the 
right hand sides of the expressions (||), while the second terms there represent 
the continuous parts of these integrals determined by the principal values of 
the Cauchy integrals with the same densities. The integral equations (|l^) arise 
immediately, if we use the integral representations mentioned just above in the 
first two algebraic constraints in each line of (|Io|). 



To explain the structure of the integral equations (13), we recall that the 
integrals there are calculated over the cut consisting of two disconnected parts: 
L = L + +L-, that each of the Cauchy type integrals in ( |l3| ) splits into the sum of 
two integrals over L + and L_ respectively and only one of these two integrals is 
a singular one. In the integrands [A]f means the jump of the function A+(£, rj, £) 
if ( = £+ G L+ or of the function A_ 77, C) if C = C- € The unknown 
vector functions <p (£, 77, t), ip (£, 77, r) and the vector functions k(r), l(r) in ( fL3| ) 
should get the suffix "+" , if their argument r = t+ 6 L+ and the suffix "— " , 
if their argument r = r_ and the corresponding suffixed vector functions 

should be identified with the fragments of the local structure of '4' defined in 
(ph. Because the parameter r in the equations (|l^) also run the entire cut L, 
i.e. it should take the values r = t + G L + as well as r = r_ G each of 
the equations (113) represents a coupled pair of vector integral equations. More 



explicitly, the first of the equations ([L3|) can be represented as a system 
— /-: <P+(C+)^C+ + — / -7 (p-(C-)dC- =k+(r + ) 



7TI 

L+ L- 



7TI 



with unknown vector functions cp + 77, t+), cp _(£, 77, r_) and the scalar kernels 

/C(r ± ,C + ) = -[A + ] f+ (k ± (r ± )-l + (C+)) 

/C(T±,C_) = -[A_] c _(k±(r±).l_(C-)). 

Similarly, the second equation in (|l3| ) leads to an equivalent system of integral 
equations for two unknown vector functions ip + (£, Vi T +)i 4" -(£) r -) : 

— /-: ip+(C+HC+ + — / -7 *P- (C-)^C- = 1+ (t+) 

1 ^1^)^ + + 1 /|^l^_( C _ )dC _i_ (r _) 

where the kernels JC±± possess the forms "almost symmetric" to that of /C±±: 

^(r ± ,C + ) = -[A; 1 ]c + (k + (C + )-l±(r ± )) 
^(r ± ,C-) = -[A: 1 ] c _(k_(C_)-l ± (r ± )). 
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Vector solutions cp +(£, 77, t+), cp_(^r?,r_) and ^ + (£, 77, t+), ip-(^,r?, r_) of 
these equations together with the corresponding monodromy data vectors k+(r+), 
k_(r_) and with expressions ( |l2"| ) for l+(r+), l_(r_) determine the solution of 
our spectral problem (^|) - (ph) by means of the quadratures 

1 r ik 1 ] 



*(Z,r ] ,w)=I+- f Liii+^ + (C + )®k + (C+)dC+ 
m J (+-w 

+-. I £=-l^t|>_(C-)®MC-)dC- 
*- 1 (£,?7,™) = i + - / -N^i+(C + )®<P + (C + )dC+ 

m J Q+ -w 

+- f -N^i_(C-)®<P-(C-)dC- 

7n 7 C- — w 

3.5 Calculation of the field components and potentials 

All components of the solution can be expressed in terms of the matrix R(£, 77) 
determined by the asymptotic expansions ]l9fl 

V = I + ur 1 Il + 0(w- 2 ), * 1 = I-w- 1 R + 0(w- 2 ) (15) 

Hence, for this matrix we have the following alternative expressions: 

R=- / [A+] c+ 1 + (C+) ® cp + (C+)dC+ + — / [A_] c _l_(C-)®«p_(C-)dC- 
m J m J 

1 /[A; 1 ] c+ ip + (c + )®k + (C + )dC+-- /[A^fc.tMC-JsMC-JdC- 

7TZ 7 7TI 7 



7TI 



The matrices U, V, W and the Ernst potentials then can be expressed as follows 



U = 2id{£l, V = 2id n TL, S = e - 2iR 3 A 

+ < ( 16 ) 

w = w (w) -4?:(OR + R t rj), $ = 2^3 5 

where the components Ra b of the 3 x 3-matrix R are numbered by the index val- 
ues A, B . . . = 3, 4, 5. The corresponding expressions for the metric components 
g a t with a, 6, . . . = 3, 4 and nonzero components of complex electromagnetic 
potential $ a are [191: 



533 = eo - i{R3 A - #3 4 ) + $3*3, 

,934 = -*(/3-/3o)+ W + i?4 4 ) + $ 3 $4, j r ; |=2/( n \ j (17) 



544 = eoe«o + i{Ri ~ Ri ) + *4$4, 




R3 5 



R 
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4 The "integral evolution equations" 

In this section we present a new way for construction of integral equations equiv- 
alent to the dynamical part of reduced Einstein equations. As in the previous 
section we concentrate our consideration on the vacuum Einstein equations and 
electrovacuum Einstein - Maxwell equations, however this construction of the 
integral equations can be realized for any of integrable reductions of Einstein 
equations mentioned in the Introduction. The structure of the derived here new 
equations called also as "integral evolution equations" differes essentially from 
the structure of singular integral equations (|l^) and their simplest regulariza- 
tions described in various forms in 1 16l, Elf. 



4.1 The "in-states" and (77,10) 

As the first step we introduce two particular values w) and \P-(r?, w) of 

the matrix function \&(£, 77, uj). In view of some analogy with the scattering 
problem we call them as "in-states" . These "in-states" are defined as 

= *(£,Vo,w), *-(»?, to) = *(£o,»7,u0- (18) 

For the hyperbolic case these are the boundary values of r],w) on the 
characteristics which pass through the reference point Po(t;o,r]o) in the orbit 
space, while for the elliptic case and are the limit values, when 77 — > 770 
or £ — * £0 respectively, of such analytical extension of the matrix function 
77, w) which arguments £ and 77 are independent complex variables instead 
of being complex conjugated to each other. The matrix functions ( |l8| ) can 
be defined also as the normalized fundamental solutions of the linear ordinary 
systems which are the restrictions of the system (^) to the corresponding char- 
acteristics 77 = ?/o and £ = £0 in the orbit space (for the hyperbolic case) or 
the restrictions of the analytically extended system (0) to the complex surfaces 
77 = 770 and £ = £0 in the complexified orbit space (for the elliptic case): 

2*(«;-o^*+ = u(f > t to ).* + r 2i(w-Ti)d v *- = v(e , »?)•*- 

*+(£o,w) =1 \ *-(»», to) =1 

Therefore, the function w) and its inverse are holomorphic outside L+, 

while *f?-(r],w) and its inverse are holomorphic outside L_ and each of these 
matrix functions posses only two brachpoints on the spectral plane, and these 
branchpoints coincide with the endpoints of the cuts L + and L_ respectively. 
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On these cuts i b± possess the structures similar to ([)]): 

L+ : #+(£, w) = A; 1 ^ +(e, w) ® h+(w) + M +(£, «;) 

= a+i+h ® cp +(e,w) + n 0+ (^w) 

L_ : *_(t7,«j) = AI^o-^.u;) ®k_(io) + M _(?7, w) 
^Z}{r] : w) = A_1_(V) <g> cp _(ry,w) + N _ (?7, «;) 



(19) 



where A + and A_, k.±(w) and l±(w) are the same as in (^); each of the row 
and column vectors tp o± , <P o± an d the matrices M ± , N ± ) are regular (holo- 
morphic) on the cuts L + or L_ respectively their suffices, and the following 
algebraic constraints on the fragments of these local structures similar to ( |l0| ) 
are satisfied on the cuts L±: 

L+ : k+ • N 0+ = 0, M 0+ • 1+ = 0, cp + ■ M 0+ =0, N 0+ • t|> + = 

L_ : k_-N _=0, M _-1_=0, cp o _.M o _=0, N„_-i|>o- = 



4.2 "Scattering" matrices x±(£> an d "dynamical" mon- 
odromy data m + (rj,w) and m_ (£,«;) 

For the next step of our construction we introduce the "dressing" or "scattering" 
matrices X±{^i 7 1t w ) presenting tj, w) in two alternative forms 

™) = x-(£)*?,w) •*-(»?> w) 

The basic idea of such alternative representation is based on the fact that the 
monodromy properties of r/, w) on the spectral plane w are conserved dur- 
ing the evolution of the fields prescribed by the reduced Einstein equations. 
Therefore, the monodromy properties of r], w) should be the same as for 
'$'+(£, w) near the cut L + and the same as for ^-(77, w) near the cut L_. This 
allows to conjecture that the matrix function x + (£,>Vi w ) should be regular on 
the cut L_|_ and X-(^Vi w ) should be regular on the cut L_. And, we have 
all what is necessary to check this conjecture. Namely, we know the analytical 
properties of \& and *!?± on the spectral plane and therefore, the structures of 
X±(£, 77, if) on the spectral plane w can be described in details. In particular, 
it is easy to see that each of the matrices X±(£,-i r ]i w ) is holomorphic function 
everywhere outside the cut L = L + + L- . Using the expressions 

X+(£,J7,«/) = *(£,77,io)-*+ 1 (£,?7o,'u0, X-(f>»7, ™) = *(C.'7>'^)-'*- 1 (^o ) ?7,io) 
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we get the local structures of x ± and X ± 011 L+ in the forms 

X+ = (k+(w) • l+0))\p + (£,?7,w) (g> cp 0+ (C,w) + M + (£,ri,w) -N 0+ (£,w) 
x; 1 = (k + (w) -l + (w))^ +0 (£,w) ® (p +(£,v,w) +M 0+ (£,w) -N+^r),™) 
X- = +(£,r),w) ®k+{w) +M+(Z,r),w)) 

XI 1 = *_(»»> «0 • (A+l+H (g><p+(£,?7,™) + N+(£, »?,«;)) 

and on L_ in the similar forms: 

X+ = (A^tp-^^^Sk-H + M-^^ic)) 

X + 1 = *+(£, to) ■ (A-l- (w) <g> cp _ (£, ?7, «;) + N_ (£, 7?, w)) 

X - = (k- (w) ■ l-(w)) o|> _ (£, ?7, w) ® <P o- (r?, w) + M_(f , 77, to) • N _ (77, to) 

Xl 1 = (k_(w) • I_(to))t|)o-(j?jt«) ® cp_(^,ry,w) + M -(?7, u>) • N_(£,T7,to) 

As it was expected, the functions x + , X+ 1 and X-i X- 1 turn out to be regular 
on L + and L_ respectively. These expressions show also that x+ and its inverse 
possess the jumps on £_ and x - and its inverse possess the jumps on L + . These 
jumps are degenerate matrices which rank is equal to 1 and which therefore, can 
be represented as the products of column and row - vectors 

[X+]i_ = [A: 1 H-(e,r?,r_)®m_(e,r_) 
[X-]l+ = [A+ 1 H+(£,rj,T + ) ® m+(?7,r + ) 

(21) 

[X + }i- = [A-]p_(£,T_)®cp_(£,r7,T_) 
[X- 1 ]l+ = [A + ]p + (?y,T + ) (8) cp + (^,77,r+) 

where t + G L + and r_ G L_. Here it is important that the vector-functions 
m_(tt, w), p_(u,w) and m+ (77,111), p + (i],w) are dependent (besides the spectral 
parameter w = r + on L + or u> = t_ on L_) upon only one of two coordinates 
£ and 77 respectively. These vector functions are determined by the expressions: 

m- (£,«;) = k_(u;)-¥+ 1 (£,u;) p_(C,«>) = ¥+(£,«;) -1_(«0 
m + (7;,7yj) = k+(iu) • ^I 1 (17,10) P+(?7,77j) = *_ (77,10) • L+(io) 

The components of these vector functions are determined completely by the "in- 
states" 10) and VP-(f?, to). It can be recalled here that the components of 
the conserved monodromy data vectors k±(io) and \±{w) in turn are determined 
by as certain fragments of the local structures ( [l9| ) of these matrix functions. 
It is useful to note here also that the vector functions (|2^) can be interpreted 
as the monodromy data for the scattering matrices X± on the spectral plane. 
To show this we consider the monodromy matrices T± which relate, similarly 
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to (pi]), the values of x± with their analytical continuations along the paths t± 
surrounding the branchpoints w = £ and w = rj in the clockwise directions: 

X-(t,V>™) — ^ X-(£,V,u>) ■ T+(t),w), X + (Lv,w) X + (£,,V,w) ■ T-(C, w) 

The described earlier local structures of x ± on the cuts L± permit to derive the 
following expressions for these new monodromy matrices 

f +M = • t + H • = i - 2 P + fa^ m + fa-) 

(p + (?7, w) • m + (r/, w)) 

T_(£,w) = *+ £,w) • T-(w) ■ * + (£,w =1-2- — r — --, 

(P-((,w) • m_(£,u>)) 

where we have used again the properties of the functions X±, that A + — ^ — A + 

and A_ — — > — A_. Unlike ( |ll| ) the matrices T± possess some evolution (i.e. a 
coordinate dependence), however this evolution do not violate their property to 

~2 

satisfy identically the relations T ± = I, which characterizes the conservation of 
the simple algebraic character of the branchpoints of x ± • (In the presence of a 
Weyl spinor field the character of the branching and therefore, the expressions 
for the monodromy matrices T± change and the identities mentioned above do 
not take place. Thus, we see that the "projective" vectors m + (r),w), m_(£, w) 
and p + (r],w), p_(£, uj) play a role of some evolving analog of the conserved 
monodromy data (or the scattering data) k±(iu), l±(w) for the solutions of the 
system (1), and hence, we call them as "dynamical monodromy data". 



4.3 Consistency of representations ( p0|) for \l>(£, ?y, w) 

Earlier it was shown, that the matrix function H?(£,r],w) is holomorphic on the 
iu-plane outside the cut L + + L_. Similarly, is holomorphic on the entire w- 
plane outside the cut L + and is holomorphic on the entire w-plane outside 
the cut L_, and the same is true also for the inverse matrices respectively. These 
properties and the local structures (|J) and ([II]) allowed us to conclude in the 
previous subsection that the function X + (£j w ) is regular on L+ and possesses 
a jump on L_, while the function X-(£i r li w ) is regular on L- and possesses 
a jump on L + . It is easy to see also that these functions should satisfy the 
conditions 

X + (£> rj,w = oo) = I, X-(€>V> w = oo) = I 

This means, that we can represent x+ an d X- as the integrals over and 
over L + respectively: 

ITT J — W ITT J C+ — w 

L_ L + 
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Substituting here the expressions ( pl| ) for the jumps [x±]i T we get the following 
representations of these matrices 



ITT J C_ — W 

X -(^V,w)=I + ^- f {^itP + (e,r y ,C+)®m + (, ? ,C+)rfC+ 
m J C,+ -w 

Similar integral representations can be derived for the inverse matrices: 

i r [a- 



(23) 



x+ 1 ^,v,w)=i + ^- f -Lip-K,C-)®«P-(e^,C-)rfC- 

tir J £_ — w 

L+ 

From the local representations (|l9|) for t + 6 L + and t_ G L_ we get: 

i+ = [*+]r+ = [A; 1 ] T+ il) +(e,T + )®k + (r + ), 
[^V = [A+] t+ 1+(t + )®(Po+(£,t + ), 

L_ : [*_] T _ = [Al^o-^T-lgk-fr-), 
I*! 1 ]^ = [A_] r _l_(r_)®cp _(r;,r_). 



(24) 



(25) 



Now we return to (|2(j) and consider a condition of consistency of these al- 
ternative expressions for This means that we have to satisfy the condition 
which we present in two equivalent forms 

X + (€,ri,™) ■ *+(£,*") =X-(f>>7>u>) • *-(»?, k>) 

(26) 

-X+ (£,»/>«>) = *_(»?, u>) -X- (f,»7.u>) 

These conditions can be reduced considerably, because both sides of each of them 
are analytical (holomorphic) functions of the spectral parameter w everywhere 
on the spectral plane outside the composed cut L = L + + L_ and both sides 
take the same value at w = oo which is a unit matrix. Therefore, the conditions 
( p6| ) are equivalent to the condition that the jumps of the left and right hand 
sides of each of the equations ((2^) on L + as well as on i_ should coincide: 

L+- X + {t+) ■ [*+]r + = [X-]r + • *-(r+) 

[*; 1 ]r + -x; 1 (r+) = *: 1 (r + ).[x: 1 ] T+ 

(27) 

L_: x - (r- = [x+]r_-*+(r-) 

[*: 1 k-x: 1 (r-) = *; 1 (r_)-[x; 1 k 
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4.4 Coupled systems of "integral evolution equations" 



Using of the expressions (|23|)- (g5|) in (|27|) leads to the following coupled pairs 
of the linear integral equations for different fragments of the local structures 
(|). One of them is a system for the vector functions ip ±(£ , 77, w) 



(28) 



^>-(^V,r-) - J 5_(^,J7,r_,C+)i|)+(f.»y,C+)dC+ =^o-(»7,t-) 

where r + , € L + and r_, G L_, and the scalar kernels 5+ and <S_ are: 

5 + (f,» 7 ,r + ,C-) = i^-^(m_K s C-)-tI»o + «,r + )) 
t_,C+) = il^+Jit. (104.(17, C+) -tDo-^r-)) 

Another (alternative) system of equations which possess the vector functions 
<P ±(£j ?7> f ) as unknown variables takes the form 



<P+(£=»7,T+) - J S+(€,TI,T + ,{-)<P -(t'Vit-) <K- = <PO+(f,T+) 



L - 



<P-(€,V,t-) - J S-(£,r),T-,(+)<P+(€,V,(+)dC+ = (Po-(t),t-) 



(29) 



where the scalar kernels 54- and S- are: 

S+(tv,r + ,C-) = - 7^- (<P 0+ (£, r+) • p_ (£, C-)) 
S_(C ) 7 ? ,r_,C+) = (<p -far_) •p + (r?,C+)) 

The structures of the derived equations ( p8|) or (|29| ) may seem to be simple 
enough, but another form of these equations can be found useful as well. 

4.5 Decoupled "integral evolution equations" 



A substitution of tp _ from the second of the equations (|28|) into the first and 
substitution of ib + from the first of the equations (p8|) into the second leads to 
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a pair of decoupled equations 

^ +(r+) - J F+{t+, C+)ip +(C+) d(+ = f+(r+) 

L+ r (30) 
t|> - (r_) - / T- (t_ , C- - (C- ) dC- - f- (r- ) 

i_ 

where the dependence of the kernels, the right hand sides and the unknown 
functions on £ and rj was omitted for brevity. The kernels and right hand sides 
of (pp|) possess more complicate than in (p8|) and d2^) structures: 



F+(t+X+) = j S+{r + ,X-)S-(x-,C+) dx-, 
f+( r +) = A>o+(t+) + / S+(T + ,x-)A>o-(X-)dx- 



f_(r_) =tp _(T_) + y <S_(r_,x+)ipo+(x+)dx+ 

Similarly, we arrive at decoupled equations for the vector functions <p ± : 
V+(t+) - J £ + (T + ,C+)cp + (C+)dC+ =S+(r+) 

cp-(r-) - j a_(r_,C-)«P-(C-)dC- =B-(r-) 

which kernels and right hand sides possess the expressions 
0+(t+,(+) = J S+(T+,X-)§-(X->C+)dX-, 

G-(r-,(-) = J l S_(r_, X+ )5 + (x + ,C-)rfx+, 



(31) 



+ ( r +) = ( Po+(t+) + / S+(t+,X-)<P o-(x-) dx- 



L- 



g_(r_) = cpo-(r-) + y 5_(r_,x+)<Po+(x+)dX+ 

For construction of any solution of reduced Einstein equations it is sufficient 
to solve only one of the systems of coupled vector integral equations (28) or (E9h 
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or their decoupled vector forms (^) or (|3l|). Besides that, all of the equations 
derived above, being vector equations, have scalar kernels and hence, they de- 
couple also into separate equations for each of the vector component. In terms 
of their solutions all components of the solution of Einstein equations can be 
determined in quadratures. 



4.6 Calculation of solution components 



In section 3.4 we recall the expressions for the metric components g a b (a,b, . . . = 
3,4), the nonzero components of a complex electromagnetic potential <i> a as 
well as the matrices U, V, W and the Ernst potentials for any solution of 
reduced Einstein or Einstein - Maxwell equations in terms of the components 
Ra b {A,B, . . .) of the matrix R defined by the asymptotic expansion ( [lEl) of 
for w — ► oo (sec ([lit), (|l7| ) and fl9fl ). Here we present the expressions for 
these components of solution in terms of the matrices R±(£, rj) defined by the 
asymptotic expansions of the scattering matrices 

X± = I + io- x R ± + 0(w~ 2 ), xZ 1 = I w _1 R± + 0(w~ 2 ) 

In particular, for the matrices U, V and W the following expressions can be 
derived easily from asymptotic considerations: 

U(£, v) = U(e, ryo) + 20 £ R+, W(£, v , w) = W(£, Vo , w) - 4i(HR+ + R^n) 

V(f, v) = V(£ , V) + 2i9„R_, = W(6, V, w) - 4*(ftR_ + R f _tt) 

For calculation of other components of the solution one can use the expressions 
( P~7|) where the matrix R defined in ( |l5|) is expressed in terms of the matrices 
R±. This last expression can be presented in two alternative forms: 

R«,»7)=R+K,t ? )+R_K > 7»)=R_K i » 7 )+R + (&,»7). 

To complete the present construction we present the alternative expressions for 
R± which follow from the asymptotic expansions of the integral representations 
(g) and (pi: 



R+ 



R 



--. / [AlV C-)®m_(e,C-) *C- 

m J 

- f [A-] c _p_(f > C-)®«P-K,»7,C-)dC- 
m J 

-- f [A; 1 ] c+ tP+(C,^C+)«'m + (ry,C + )dC+ 

— / I A +k+ p+(^C+) ®«p+(^»?>C+) <^C+ 
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These expressions allow to calculate all components of solutions also in the 
alternative forms. For example, for the Ernst potentials we have 

£(Z,v)=£(t,rio) + ^ J [AlVtei-tMC-)) (m_(C_) • e 2 )dC- = 

L_ 

= + - J [A; 1 ] c+ (e 1 -iMC+)) (m+(C+)-e 2 )dC- 

mv) = mm)-^ J [AlV(ei-tlMC-)) (m_(C_)-e 3 )dC- - 
= mo,v) -\ J [A; 1 ] <+ (e 1 .ip+(C+)) (m + (C+)-e 3 )dC- 

where ei = {1, 0, 0}, e 2 = {0, 1, 0} and e 3 = {0, 0, 1}. 

5 Concluding remarks 

In this paper the two-dimensional space-time symmetry reductions of vacuum 
Einstein equations and electrovacuum Einstein - Maxwell equations have been 
presented in some new linear (quasi-Fredholm) integral equation forms. Similar 
equations can be derived for all other known integrable reductions of Einstein 
equations using the same method without any its essential modifications. 

The alternative representations (|20| ) of the solution & of associated linear 
system in terms of "in-states" and the scattering matrices x± used here may 
be considered as some analogue of the well known dressing methods, developed 
for solution of various completely integrable systems (see [^6|, |27], [| and the 
references there). In some points the present construction can remind also 
Krichever's construction pi| of the "analogue of d'Alembert formula" for the 
Sine-Gordon equation, as well as a construction of the homogeneous Hilbert 
problem and the matrix linear integral equation form of vacuum Ernst equation 
presented by Hauser and Ernst ||. In particular, in || the matrices identical to 
\&± had been used as important elements of the developed method. However, 
unlike the mentioned above constructions, closely related with formulations of 
various matrix Riemann or Riemann - Hilbert problems, the present analysis 
is based on a more detail consideration of some features of the structure of 
reduced Einstein equations which allow to reduce the problem to much more 
simple, scalar quasi-Fredholm integral equations. 

In comparison with the previously derived (in the framework of the same, 
monodromy transform approach) singular integral equation form of reduced 
Einstein equations or their regularizations, which were expressed in terms of 
conserved monodromy data |fl8| , [HJ, the new integral equations are designed 
mainly for consideration of initial and boundary value problems for the reduced 
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Einstein equatons. Though the "old" equations had already provided us with 
a principle scheme for solution of the characteristic initial value problems and 
the Cauchy problems in the hyperbolic cases and some boundary problems in 



the elliptic cases of integrable reductions of Einstein equations |19, E9|, the 



new integral equations are obviously better adapted at least for solution of 
characteristic initial value problems. The scalar kernels and the right hand sides 
of the new equations carry more explicit information about characteristic initial 
data, because they are expressed in terms of "dynamical" monodromy data 
evolving along the characteristics and being determined completely in terms of 
a given characteristic initial data for the fields. 

It is interesting to note, that the applicability of the presented here the "in- 
tegral evolution equation" form of reduced Einstein equations is restricted from 
the beginning by the condition that the analytical (in terms of the coordinates £ 
and 77) local (i.e. near the point where the boundary characteristics intersect) so- 
lutions are considered only. This restriction leaves out of our consideration some 
kinds of characteristic initial value problems which correspond to physically in- 
teresting enough situations such as, for example, a collision of plane gravitational 
waves propagating with distinct wavefronts on the Minkowski background. In 
this case, the regularity of solutions which take place in some globally defined 
coordinates is not compatible with their analyticity in terms of the "geometri- 
cally defined coordinates" £ and rj near the point of the wavefronts collision J3(| . 
Therefore, to consider a collision of such waves we need to refine essentially on 
our methods. Fortunately the monodromy transform approach, which farther 
development was presented in this paper, admits an appropriate generalization. 
This generalization was found very recently in the author's collaborated paper 
with J.B.Griffiths |3l| where the generalized system of linear "integral evolu- 
tion equations" was derived and a method for direct solving of the characteristic 
initial value problem for given characteristic initial data for colliding plane grav- 
itational or gravitational and electromagnetic waves propagating with distinct 
wavefronts on the Minkowski background was presented. 

It is necessary to note here also, that in the cases of a Cauchy problem for 
hyperbolic reductions or boundary problems for elliptic reductions of Einstein 
equations, the interrelations between the initial or boundary data for the fields 
and the functional parameters in the kernels and coefficients of the integral 
equations derived in the present paper turn out to be more complicate than 
in the characteristic case. Therefore, for solutions of these problems there is 
no yet a similar more simple way than suggested in | fl9[ p9| general scheme 
for consideration of such problems using the linear singular integral equations 
whose construction is based on the conserved monodormy data. 
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